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Abstract 



The study of wavepacket revivals is extended to the case of Hamiltoni- 

ans which are made time-dependent through the adiabatic cycling of some 

parameters. It is shown that the quantal geometric phase {Berry's phase) 

causes the revived packet to be displaced along the classical trajectory, by 

an amount equal to the classical geometric phase (Hannay's angle), in one 

degree of freedom. A physical example illustrating this effect in three degrees 

of freedom is mentioned. 

PACS numbers: 03.65.Bz, 31.50.+W 
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Ordinarily, a quantal wavepacket following a classical trajectory spreads irreversibly, until 
it completely loses its integrity as a localized packet. However, when the underlying classical 
trajectory is stable and periodic, the dispersed wavepacket eventually "puts itself together 
again" and for a time continues, localized, along the classical trajectory. The argument for 
this comes from an examination of the dynamical phases exp —iE n t/h acquired by the energy 
eigenstates whose superposition makes up the wavepacket. |[l|-[3[ Such wavepacket revivals - 
which have been observed experimentally in Rydberg atoms [f§] as well as one-atom masers f| 
- are intriguing because they represent the seemingly spontaneous resurrection of classical 
behaviour in a quantal system. The study of wavepacket revivals has so far been restricted to 
systems where the Hamiltonian is time- independent. The present paper extends the theory 
to the case of slowly time-dependent Hamiltonians, where, as will be shown, Berry's phase 
adds a new and interesting twist to the picture. 

Specifically, in this paper I consider wavepacket evolution under an adiabatically cycled 
Hamiltonian, that is, a parameter-dependent Hamiltonian if (R), where the parameter R is 
made to slowly trace out a closed loop in parameter space. I will restrict myself to systems of 
one degree of freedom, although parameter space is multi-dimensional. Thus the underlying 
classical motion (assumed bounded) is periodic for fixed R. When R is made to change 
slowly with time, the classical motion is nearly periodic: a trajectory goes round and round 
a slowly changing closed curve in the two-dimensional phase space. At any instant, this curve 
is an energy shell — a level surface — of the instantaneous Hamiltonian, and is determined 
by the requirement that the action §pdq over one "period" is an invariant (see e.g. Ref. 
for details). The analysis that follows shows that, as in the time-independent case, there 
will be a revival of the wavepacket: the standard argument for revivals is easily extended 
to the case of a slowly time-dependent Hamiltonian, only the dynamical phases now have 
the form exp — (i/h) J E n (t) dt. However, this is only part of the story. As demonstrated by 
Berry [[7], an energy eigenstate evolving under an adiabatically cycled Hamiltonian acquires 
not only a dynamical phase, but also a geometric phase, determined by the loop traced 
out in R-space. A proper analysis of wavepacket evolution in this situation must take this 
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phase into account. I show that, if the time at which the parameter returns to its initial 
value is chosen to coincide with the revival time, then the net effect of Berry's phase is to 
cause a displacement of the location at which the revived wavepacket appears. Moreover, 
the amount of this displacement is given by Hannay's angle ||, the classical analogue of 
Berry's phase. As the analysis will show, this effect is generic in one degree of freedom; at 
the end of the paper I briefly mention a physical example in three degrees of freedom which 
illustrates the same effect. 

To demonstrate what has been stated above, it is useful to first review a result concerning 
wavepacket evolution under a time- independent Hamiltonian H; this result is summarized 
by Eq.| below. Imagine a quantal wavepacket ipo localized around some point (qo,po) in 
phase space. (That is, the configuration-space representation is localized around q , the 
momentum-space representation around Pq. More formally, one could represent the quantal 
state directly in phase space, using a Wigner transform.) Expressing ip as a superposition 
of energy eigenstates \u n ), the non- vanishing coefficients a n are confined to some finite range 
of values of n, centered, say, around n. Assume the semiclassical regime: n ^> 1. Under time 
evolution, the coefficients in this superposition acquire phases: ip(t) — J2n exp(ia n )a n |w„), 
where a n = —E n t/h. It is the process of "de-phasing" — the fact that the a n 's grow at 
different rates — which determines how the wavepacket moves and spreads. The Correspon- 
dence Principle, however, dictates that in the semiclassical limit (n — > oo), the packet is 
propelled along a classical trajectory. Let us consider this more closely. 

A Taylor expansion of E n around n yields: 

an = ~ An^E' n t - ^-(An) 2 ^ , (1) 

where An = n—n, and E' n = dE n /dn, (and similarly for the second derivative E£, and higher 
derivatives E^' ||.) Using 1/n as an ordering parameter, E^' is 0(l/n s ), and h is 0(1 /n) 
(in comparison with the relevant "macroscopic" quantities E n and I n = (1/2-k) § pdq, with 
the latter calculated for a classical orbit of energy E n ). The first term on the right side of 
Eq.|l], —E n t/Ti (which scales as n), contributes an external phase to the wavepacket. The 
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second term (~ n°) describes the linear (in An) de-phasing between the components \u n ) of 
tpo, the third term (~ 1/n) the quadratic de-phasing, and so forth. 

In the semiclassical limit, in which the wavepacket follows a classical trajectory, the first 
two terms dominate: 

ip(t) = exp(—iE n t/h) exp(— iAnE' fl t/h)a n \u n ) 

n 

= e^^^^nlUn). (2) 

n 

Thus, disregarding the external phase a n , the phases 

/3 n = -AnE' n t/h (3) 

specify a wavepacket that evolves along a classical trajectory. (Note in particular that this 
evolution is periodic.) This result is valid for times of order unity: t ~ n° (by which is meant 
that the ratio of t to the relevant macroscopic time scale — the period of an orbit of energy 
E n — is held fixed as n — > oo). At longer times, the quadratic and higher-order de-phasing 
terms become appreciable, and are responsible for the spreading of the packet. Now, WKB 
theory gives us 

U „ ™ (4) 
ft " 81 ' W 

where dH/dl is the derivative of the classical Hamiltonian with respect to the classical 
action, evaluated at energy E n . The correction to Eq.|] is 0(l/n 2 ) (the 0(1 /n) correction 
is identically zero). In terms of action-angle variables (0,1) [fj, dH/dl is simply the rate 
of change of the angle variable 9 along a classical trajectory. Thus, neglecting the 0(1/ n 2 ) 
correction to Eq||, we get 

(3 n = -AnAO, (5) 

where AO = (dH/dI)t is the change in 9 associated with classical evolution for time t. We 
therefore conclude that, if we modify a semiclassical wavepacket by tacking on linear 
de-phasing factors exp(— iAnAO) to the terms in the superposition, then the effect is to 



shift the location of the wavepacket along a classical trajectory, by a change in the angle 
variable equal to the amount of (linear) de-phasing, A9. Symbolically, 

{ Qn _> e - iAnA0 a n } (e,I)^(e + A9,I). (6) 

Let us now move on to adiabatically driven systems. Consider a parameter-dependent 
Hamiltonian H(R), and imagine again an initial wavepacket ipo localized at (<7o>Po)- Now, 
however, let the packet evolve under the time- dependent Hamiltonian obtained by making 
R slowly trace out a closed loop T in parameter space. Denote the initial (and final) point 
on this loop by Ro, and the initial and final times by t = and t = Tp. 

Assuming dH/dt slow enough that the quantal adiabatic theorem holds, the wavefunction 
at time t is 

m = £e^a n MRW)>, (7) 

n 

where the n 's are real, the coefficients a n are time- independent, and the |w n (R))'s are the 
eigenstates of if (R). Thus, when R(t) returns to R at time T r , we have a state identical 
to ipQ, except that the expansion coefficients have acquired phases: a n — > e i?in(Tr )a n . 

Before the discovery of Berry's phase 0, one might have guessed that the phases 4> n (T r ) 
are given by 

0n(T r ) = -\ f Tr dtE n (t), (8) 

h Jo 

where E n (t) is the nth eigenstate of H(R,(t)). For the time being, assume Eq.|^ is correct. 
Expanding as in Eq.|I| above, we get 

MTr) = -?SSk) - AnhSll - (A n )»^ - ■ ■ , (9) 

where 

b s (T r ) = ( TV dtE^\t). (10) 

J 

The sizes of the terms in Eq.^ are determined by a competition between the largeness of 
and the smallness of 1/n. For instance, suppose we take n — > oo, while letting Tp scale 
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as unity: Tr ~ n° (in the sense denned earlier). This is the limit in which classical behavior 
persists over the entire time of observation. On the other hand, in this limit b s (Tr) ~ l/n s , 
so the quadratic and higher-order de-phasing terms in Eq.^ are negligible. Thus, as in the 
time-independent case, we associate linear de-phasing with classical propagation, and the 
higher-order terms with spreading. 

The regime of interest in this paper is that for which Tr is large enough that the first 
three terms of Eq.|9] are significant, while terms of order (An) 3 and higher are not. (Formally: 
n — > oo, Tr ~ n.) In this regime the spreading of the wavepacket becomes important, and 
wavepacket revivals appear. 

As mentioned, a revival occurs when a packet — after having lost its localized structure 
- reassembles itself, continues for a while along a classical trajectory, then again spreads 
and dissolves. (Such revivals repeat themselves regularly, in a pattern of "quantum beats" 
Additionally, there are fractional revivals f|], but these will not be discussed here.) To 
show that revivals appear in adiabatically cycled systems, suppose we choose T r so that 
& 2 (Tr)/2/i = 2tt. Then at time Tr the effect of the quadratic de-phasing term in Eq.[5] will 
be null, and the wavefunction will reflect only the external phase exp(i<j) n ) and the linear 
de-phasing term 

(We have discarded a correction that scales like 1/fi.) Here dH/dl at time t is evaluated at 
the energy shell of if (R(t)) which corresponds to the eigenstate n, in other words at constant 
J, even as H slowly changes. According to Eq.|| Eq.[Ll] describes a wavepacket which has 
shifted from its original position in phase space, (qo,Po), along a classical trajectory of 
H(Hq), by a change in the angle variable given by 

A9 = l * w (t). (12) 

Eq.O seems to place the revived wave packet at the point in phase space which the 



classical trajectory would have reached after time Tr. That is, since 9 = dH/dl in the 



time-independent case, and since in the adiabatic limit the action / is an invariant ||10|| , it 
appears at first glance evident that the classical change in 9 after time Tr should be given 
by Eq.[12], with dH/dl evaluated at constant /. This picture, however, is erroneous: as 
demonstrated by Hannay ||, the classical trajectory experiences an additional shift in 9: 

A9 = | Q Tr dt^it) + A9 H . (13) 

This "extra" shift A9h, Hannay 's angle, is geometric: it is determined by the loop V in 
parameter space. I will now show that the source of the discrepancy — the reason Hannay's 
angle did not appear in Eq.|l2| — is the neglect of Berry's phase in Eq|| 

In Ref. , Berry showed that the correct phase acquired by an adiabatically cycled eigen- 
state consists of both the dynamical phase of Eq|| and a geometric phase 7„ (determined 
by the loop T): 

MTr) = -t I r dtE n (t) + 7 n- (14) 
n Jo 

Calculating the effect of Berry's phase on the evolution of a wavepacket launched from 
(qo,Po) is simple. First, Eq.^| gathers the extra terms 

^ + A "t + >) 2 ^ + -- < 15 > 

Now, 7^ ~ n (see e.g. Eq.[30] of Ref. [|llj), so d Sr y n /dn s ~ n 1_s . Thus, in the limit n — > oo 
only the first two terms in Eq.|15| survive. Let us again choose Tp so that &2(^~r)/2fr = 2tt. 
In the notation of Eq.|2], the wavepacket at time Tr is then given by: 

V(T r ) = e^^e^K), (16) 



where now 



1 r T T 

OLn = - 



-( T dtE n (t) + ln (17) 
n Jo 

Pn = -An[[ TT dt^(t) - ^1. (18) 
L Jo ol dn J 

This describes a wavepacket whose position in phase space is specified by both the "dynam- 
ical" shift of Eq.|l2l, and an extra shift —d^njdn in the angle variable: 



To complete the connection with the classical result, Eq.[l3], we invoke the central result of 
Ref. [PH , which states that Hannay's angle A6h and Berry's phase j n are related semiclas- 
sically by: 

M H = (20, 

Thus, the revival described by Eqs.[l6| - [18] does indeed appear where a proper classical anal- 
ysis suggests it ought to. Specifically, the dynamical quantal phases / E n dt are responsible 
for the "dynamical" shift in the angle variable (EqJfJD, while the geometric quantal phases 
7 n further boost the packet by an amount equal to the classical geometric shift A6h- 

(A word about the sizes of terms in Eq.[19|. Although the first term on the right side 
scales like n, whereas the second term is order unity, the two nevertheless have a comparable 
effect on the location of the wavepacket, since for this purpose their values are only relevant 
modulo 2tt. On the other hand, the leading-order correction to the first term scales like l/n 
- since corrections to Eq.|] are 0(l/n 2 ) — thus its effect is genuinely small.) 

Eq.|TJ5] embodies the central result of this paper; combined with Eq.|2D|, it reveals the effect 
of Berry's phase on wavepacket revivals in adiabatically cycled systems. In addition to this 
result, the formalism used in this paper provides a simple interpretation of the semi classical 
relationship between Berry's phase and Hannay's angle (Eq.pOf). Namely, if we modify a 
wavepacket i[) = J2n a n |w n (Ro)) by tacking on Berry's phases exp to the terms in the 
superposition, then the effect is to shift the packet along a classical trajectory of if(Ro), 
by a change in the angle variable equal to Hannay's angle, A9h- This follows directly from 
Eqs.|] and |2(], without any need for a discussion of revivals. (Alternatively, one could stand 



the argument on its head and derive Eq.|2(] by considering the semiclassical evolution of a 
wavepacket under H(R(t)) and invoking Eq.|6|.) This "wavepacket interpretation" of Eq.|20| 
is similar in spirit to arguments presented by Hannay ; its novelty resides in that it offers 
an easy visualization of the relationship between the two geometric quantities, Berry's phase 
(quantal) and Hannay's angle (classical). 

8 



It has been assumed in the preceding analysis that the value of Tr is chosen so that a 
revival occurs just as R(t) concludes its circuit in parameter space. Suppose now that we are 
less restrictive with the value of Tr, but still within the regime of wavepacket revivals (i.e. 
&2(Tr)/2ft is order unity, but not necessarily 2n). A straightforward calculation reveals that 
a revival then occurs (along the classical trajectory) at time Tr satisfying 5 2 (Tr)/2/i = 2tt, 
even though the Hamiltonian at t = Tr is different from that at t = 0. Thus, revivals are 
as generic in adiabatically driven systems as in time-independent ones. However, only when 
the revival time Tr coincides with the cycling time Tp is it meaningful to discuss the effect 
of Berry's phase on the revived wavepacket, since Berry's phase (as well as Hannay's angle) 
is well-defined only for closed circuits in parameter space. 

The effect of Berry's phase on the revival location of wavepackets has been illustrated 
here in one degree of freedom. It would of course be very desirable to find a physical (three- 
dimensional) system which exhibits this effect. One candidate, suggested in conversation 
by M. Nauenberg, is a Rydberg atom in a weak magnetic field B; here the classical motion 
is characterized by the precession of the Kepler orbit at the Larmor frequency, and with a 
proper choice of |B| a wavepacket launched upon the orbit will experience a revival |3|]. Now 
suppose that, upon launching the packet, we let the direction of B adiabatically trace out a 
loop which encloses a solid angle Q in B-space, so that at the revival time the field is back 
to its initial orientation. Then as a result of Berry's phase, the revived wavepacket appears 
at a location which differs from where it would have appeared had B remained constant. 
The difference is simply a rotation by Q Jl2] around the initial direction of B. |13| This is 
essentially an atomic version of the Foucault pendulum, and if experimentally observable 
|14| , would constitute a vivid demonstration of the geometric phase effects illustrated in the 



present paper. 
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